(12) 



UK Patent Application o 9 ,GB m>2 302 777 n 3 ,A 



(43) Date of A Publication 29.01.1997 



(21) Application No 9513086.0 

(22) Date of Filing 27.06.1995 



(71) Applicant(s) 

Motorola Israel Limited 

(Incorporated in Israel) 

16 Kremenetski Street, Tel Aviv 67899. Israel 

(72) Inventor(s) 

Eliezer Fogel 

(74) Agent and/or Address for Service 
Hugh Dunlop 

Motorola Limited, European Intellectual PropertY 
Operation, Midpoint, Alencon Link. BASINGSTOKE, 
Hampshire, RG21 7 PL United Kingdom 



(51) INTCL 6 

H04L 25/03 , H04B 7/005 

(52) UK CL (Edition O) 

H4P PRE 



(56) Documents Cited 
EP 0096802 A2 



US 3864632 A 



(58) Field of Search 

UK CL (Edition N ) H4P PAN POP PEF PRE PRR 
INT CL 6 H04B 1/10 7/005 7/212 . H04L 25/03 
Online: WPI 



(54) Recovering symbols of a digitally modulated radio signal 



(57) A method of recovering a finite duration of N symbols of a digitally-modulated signal, wherein a 
received signal is sampled at a period J v generating a batch of samples, wherein K, , is chose * to 1 be an 
integer multiple of the N symbols. The K 1 time-domain samples are input to a discrete *™™ r }™^™ 1 ™ 
providing a first batch of K, frequency-domain samples. The first batch of K, 1 W™**^ 
multiplied by a set of frequency-domain K n complex gains to provide a second batch of K, ^ u «^^^" 
samples. The second batch of K, frequency-domain samples are broken into ^/N batches of N samples. The 
Kl /N batches are added, generating a single batch of N samples which is input to an 

transform (24) to obtain N time-domain samples. A zero padding option is descnbed. Application is to TDMA 
systems, with compensation for inter-symbol interference. 
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filed was informal and the print reproduced here is taken from a later filed formal copy. 
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METHOD OF RECOVERING SYMBOLS OF A DIGITALLY 
MODULATED RADIO SIGNAL 

Field of th e Invention 

This invention relates to a method of recovering symbols of a digitally 
modulated radio signal. 

i 

Background of the I nvention 



In digital communications systems using Time Division Multiple 
Access (TDMA) schemes, to apportion the Available specfrum to the 
6ommunication users, digital receiver^ are typically designed for operation 
in the time-domain to recover the digitally transmitted information. As the 

15 technology used in such digital* communication's systems has increased in 
complexity, the requirement on the computational capabilities, of 'digi^d 
radios has increased. Consequently, there is an increasing need for more 
optimal use of the digital signal processing resource within digital radios. 

In a TDMA system data recovery isgsnerally carried out in the time- 

.20 domain for a whole slot of data using, for example, a Viterbi decoder. The 
traditional approach to the implementation of an optimal digital receiver 
uses a matched filter followed by a (infinite optimal or finite sub-optimal) 
tap delay equaliser operating in the time-domain. The ; matched filter at the 
receiver corresponds to the transmit filter of the digital transmitter. In 

25 digital systems the transmitted data symbols are convolved with the 

transmit ^filter impulse response (git)) to provide the transmitted signal, s(t) 
= g(t)* data. This transmitted signal needs to be demodulated and the 

j original data symbols recovered at the digital receiver. 

A time-domain implementation typically uses an equaliser to 

30 compensate for any iriter-symbol interference (ISI) of the received symbols 
that has been caused within the channel. This is not necessarily an optimal 
.solution. Moreover, having a complex channel, impulse (in the sense of real 
and imaginary parts), requires a complex finite impulse response (FIR) filter 
for equalisation purposes, making the solution even more computationally 

35 - expensive. - - * . - - - 

Typically, such operations are computationally complex, requiring a 
significant amount of the signal processing resource of a digital radio. 
Additionally, in order to accommodate the use of time-domain equalisers, 
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effort is directed at achieving a Nyquist piilse : shape response for the 
transmitted and received symbols/ 

Thus there is a need for an improved method of recovering symbols of 
a digitally modulated radio signal which avoids some of the disadvantages of 
prior art arrangements. 

Summary of the Invention " ' , : 

In a first aspect of the present invention, a method of recovering a 
finite duration of N symbols of a finite duration digitally-modulated received 
radio sighalis provided: The method includes the steps of receiving a signal 
ahd sampling the signal at a sampling period T-i to generate a batch of Ki 
time-domain samples, wherein the sampling frequency is a multiple of the 
symbol frequency. The batch of Ki time-domain samples are transformed to 
the frequency-domain to provide a first batch of Ki frequency-domain 
samples: The first batch of Ki frequency-domain samples are multiplied by 
a set of frequency-domain Ki complex gainst (ai: ix=l: to Ki), to provide a 
second batch of Ki frequency-domain samples. The second batch of Ki 
frequency-domain samples is broken into Ki/N batches of N samples which 
are added to each other to generate a single batch of N samples. The single 
batch of N samples id inverse transformed to obtain N time-domain samples 
which are the desired transmitted symbols'.'-' : 

In this maimer, a batch of (Ki) digitally-modulated received time- 
domain saihples, whereinKi is chosen to be a multiple of the number of 
symbols N and the sampling frequency is a multiple of the symbol frequency 
of the received signal, are transformed into the frequency-domain, processed 
and inverse transformed back into the time-domain to recover the desired 
symbols. r 

In a second aspect of the present invention, a method of recovering a 
finite duration of N symbols of a finite duration digitally-modulated received 
' i-adio signal is provided. The method includes the steps of receiving a signal 
and sampling the signal at a sampling period Ti to generate a batch of Ki 
time-domain samples, wherein the sampling frequency is not a multiple of 
the symbol frequency. The batch of Ki time-domain, samples are 
transformed to the frequency-domain to provide a batch of Ki frequency- 
domain samples. Zeros are appended to the Ki frequency-domain batch of 
samples to provide a batch of K2 samples in the frequency-domain where K2 
> Ki. The batch of K2 frequency-domain samples are multiplied by a set of 
frequency-domain K2 complex gains (af: i=l to K2), to provide a batch of K2 
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frequency-domain samples. The batch of K2 frequency-domain samples is 
broken into K2/N : batches of N samples which are added to each other to 
generate a single batch of N samples. The,single batch of N samples is 
inverse transformed to obtain N timerdomain samples which are the desired 
5 transmitted symbols. 

In this maimer, a batch of digitally-modulated received time-domain 
samples, wherein the sampling frequency is not a multiple of the symbol 
frequency, can be recovered by transforming the batqh into the frequency- 
domain, processing the frequency-domain samples and inverse transforming 

10 the batch ; back into the time-domain tq recover -the desired received time- 
domain symbols, K2 is assume^ to bejpreajter than Kl with t}ie 1^1 

, , frequency-domain /samples are .padded by (K2 - Kl) zeros. The zero padded 
frequency? domain samp Jes correspond tq the sampling rate of the received 

, : signal; which is~a mjiltiplipative of N t tike number of symbols 

15 In 3 preferred, embodiment of the; present invention, the K2 complex 

. gains tei: i=lto.K2) K are calculated.by evaluating the receiver pulse shape 
filter response giin the frequency c (^main to c obtain the batch ,of K2 complex 
values (gi: i^l to, K2>- Further, the complex gains gi (gi:.i=l to K2) are 
normalised by breaking the i^atch of K2 complex gains gi into K2/N batches 

20 of N&amples (gij: i=l.taN^l tpjC^N) and adding the absolute squares of 
the K2/N batches to each other to generate^ single batch of N normalising 
. v values (qi: i=l to N) such that qj=.sum[ I gij 1 2 ; j =i to K2/N]. Constants ei 
may be added to each qi to alloyr the introduction of arbitrary design 
parameters into the frequency-domain receiver ^design process. Normalising 

25 (gij: j=l to N) of each j is then performed by dividing the values gij by the 

values qi to obtain a new set of gains (aij=gij/qi : i=l to N, j=l to K2/N). The 
K2/N subgroups (aij) are thei* combined back into a single group to be used 

, as the (ai) complex gains. ^ ti _ ; 

In this manner, the computational complexity of calculating the 

30 inverse transfer of the batch off frequency-domain samples is greatly 

- , simplified. - * • t . 

T^e^ecei^er pxilse,s}iape. filter response gi may be pre-stored, or 
alternatively the parameters aj may be pre-calculated and pre-stored if the 
channel characteristics are.known apriori. - . . ^ ' 

35 A preferred embodiment of the invention will now be described, by 

. way of example only, with reference to the drawings. , - 
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FIG. 1 is a block diagram of the signal processing functions of a prior 
art digital radio receiver- r ' — 

FIG. 2 is a block diagram of a frequency-domain digital radio receiver 
in accordance with a preferred embodiment of the invention. 

FIG. 3 is a flow diagram illustrating a method of recovering a finite 
duration of N symbols of a digitally-modulated received radio signal in 
accordance with the preferred embodiment of the invention. 

FIG. 4 is a flow diagram for determining the complex gain values aj, 
from the receiver pulse shape filter response, in accordance with the 
preferred embodiment of the invention. ' • , , i 

Detailed Descrintion of the Drawings • * - : 

Referring first to FIG. 1, a block diagram of the signal processing 
functions of a prior art digital radio receiver is shown. The signal processing 
functions comprise a 1 receiver matched filtfer g(-t) 10, an equaliser h(t) 12, a 
sampler 14 and a quantizer 16. - < » ■ * . 

In operation, the digital receiver receives a time-domain signal x(t) 
which is input to a matched -filter g(-t) 10. The matched filter in the receiver 
is a time reversed replica of the puis** shape filter in the transmitter. 
Implementing a matched filter 10 in the receiver ensures that the 
combination of both filters gives the desired overall pulse shape response. 
The filtered pulse is then input to an equaliser h(t) 12 to remove the channel 
effects of inter-symbol interference (ISI). In a practical time-domain 
receiver, abd given a realistic rebeived pulse, ISI is introduced in the 
hardware filtering and channel filtering. If the equaliser 12 is a linear 
equaliser the equaliser operation may be performed before the matched 
filter 10. The ISI-free received pulse is then input into a sampler 1*4 to 
sample the signal at a sampling period T, generating a batch of N time- 
domain samples. The batch of time-domain samples is ; then input into a 
quantizer 16 to obtain the desired symbols. * * 

Specifically, consider the case of a finite duration signal i.e. a time slot 
in the TDMA context, consisting of N symbols s(n): n=0...N-l and let the 
noise free received signal be: - \ - i * 

N-l 

x(t)= £s(n)g(t-nT) [1] 

n=0 
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where 

s(n) is the n-th symbol ; * _ - ' * : . - 

g(t) is the pulse shape 
T . is jthe symbol period 

then the optimal detected symbol is expressed as 



s(n) = Q 



X^gRn-iq)T-t)*x(t)J 



[2] 



-Where* is the convolution operator, 
Q[ ] is "quantization" into the, symbols alpha-bet, and q£ 1 
,10. ci are the equalizer coefficients , : ^ . , *_ * . 

The choice of the coefficients, Jci) depends on the optimisation 
criterion. q<l corresponds to a fractionally spaced equalizer (FSE) chosen so 
that l/(qT) is greater than the signal, bandwidth, as known to those skilled 
15 in the art. The number of coefficients {cil and their domain (real or complex) 
depend on the channel characteristics. ■ 

In equations [1] & [2] time (t) may be interpreted as a continuous 
. variable or as an integer ifx(t) is a sampled signal. In the present context, t 
is considered an integer. The physical implication is that x(t) is the t-th 
20 sample of. a, continuous signal, sampled at an appropriate sampling ^jate 

1/Ti. Note that T in general, for a, sampled signal, is not riecessaril^an 

integer, i.e. Ti/T is not.necessarily an integer. This has important ^ 
ramifications in the practical implementation of the, receiver, described 

- # herein. • \../> . ■ .» : • .1 . * - 

25 - The prior art for digital receiver implementation is summarised in 

f 4 equation~[2]. For practical channels, the receiver in equation [2] requires 
- - : time, phase and frequencysynchronisation, in addition to the convolution 
r and equalization operation of equation [2].. r ; v , v . s , [ 

In typical time-domain digital receivers, where the time-domain 
30. symbols are processed sequentially, the signal processing. needed to 

implement the equaliser and matched filter operations is complex and 
computationally intensive. , r - . vS 

Referring now to FIG.2, a block diagram of a fi-equency-domain 
digital radio receiver, in accordance t with a preferred embodiment of the 
* 35 invention, is. shown.- The signal prpcessing. functions comprise a discrete 
fourier transform (DFT) 20 with optional zero padding, an overlap add 
function 22, an inverse fast fourier transform (IFFT) 24 and a quantizer 26. 
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In operation, the digital frequency-domain receiver receives a time- 
domain batch of Ki samples x(t) which are input to the DFT 20. The DFT 
20 transforms the batch of Ki time-domain samples into the frequency- 
domain to provide a batch of Ki frequency-domain samples. The appending 
5 of zeros to the Ki frequency-domain samples can be implemented at this 

stage to provide K2 frequency-domain samples if the sampling frequency 1/T 
is not an integer multiple of the*symbol frequency 1/Ti, e.g. T/Ti is an 
integer. If the sampling frequency is an integer multiple of the symbol 
frequency, Ki = K2. The received samples in the frequency-domain are then 
10* input to an overlap add function 22. The overlap add function 22 multiplies 
the batch of K2 frequency-domain samples by a set of K2: complex gains (ai: 
i=l to K2), to provide a batch of K2 frequency-domain samples. The batch of 
K2 frequency-domain samples, are then broken into K2/N batches of N 
samples which are added to each other to generate a single batch of N 
15 samples. The frequency-domain signal is then input to an IFFT 24 to 

transform the single* batch of N samples back into the timer domain to obtain 
N samples. The N samples are then input into a quantizer 26 to obtain the 
"symbols" output of the received signal s(n). 

Referring now to FIG; 3, a flow; diagram is shown illustrating the 
20 operation of the receiver of FIG. 2 in greater detail. A time-domain batch of 
Ki pulses x(t) are input to an analog-to-digital converter, as in step 30, to 
provide Ki time-domain digitised samples as shown in step 32. The digital 
batch of K'l time-domain samples is input to a DFT 20, as in step 34. The 
DFT 20 transforms the batch of K.J time-domain digitised samples into the 
25 frequency-domain to provide a batch of Ki frequency-domain samples. 
Zieros are appended to the Ki frequency-domain samples to provide K2 
frequency-domain samples if the sampling frequency is not an integer 
multiple of the symbol period^ as in step 36. If the sampling frequency is an 
integer multiplt/of the symbol frequency, Ki = K2. The K2 frequency- 
30 domain samples are multiplied, :in step 38, by a set of K2 complex gains (ai: 
i=l to K2) that represent the receiver pulse shape filter response in the 
" : frequency-domain, as shown in step 40, to provide a batch of K2 frequency- 
domain samples. The batch ofXg frequency-domain samples is then broken 
into K2/N batches of N samples, as shown in steps 42 to .46, and these 
35 batches are added to each other (step 48), to generate a single batch of N 

samples, as in step 50. The process of adding the K2/N batches of N samples 
to each other is equivalent to overlapping the K2 samples to. generate a 
single batch of N samples. The value of N dictates the number of additional 
samples to be generated between step 44 and step 46 of FIG. 3. The N 
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frequency-domain samples are then input to an EFFT, as in step 52, to 
transform the single batch of N samples back into the time-domain. The N 
time-domain samples are then.input.into a quantizer 26, as in step 54, to 
obtain the "symbols" output of the received signal s(n). . 
5 A frequency domain implementation of ithe symbol recovery is 

: especially suited for TDMA systems as an entif e.slot is analysed in "batch" 
processing compared to a time domain implementation where conventional 
' sequential processing is performed j Performing the batch processing in the 
frequency-domain removes the constraint on the pulse shape to be ISI free. 
10 Advantageously the pulse shape can be designed to accommodate other 
; ■ system xohsideratione such. as inter-cha^el interference (ICD- In addition 
to removing the piifce shape constraint, channel .time invariant frequency 
distortion^ removed as the frequency-domam-receivertis optimally 
. < recovering the symbols from^tha received pulse shape. The radio frequency 
15 . (RF) channel and non-symmetric intermediate frequency (IF) filtering are 
7 then absorbed trivially' into the receiver design. Filtering e.g. to eliminate 
narrow band disturbances from ^jctejrnal or internal sources such as carrier 
feedthrough, can^be applied withlittle effect on performance and 
insignificant computational -load. - >* * , . 
... 20 Referring now to FIGi 4, a flow diagram for calculating the complex 
, ■ gains ai.from.the received pulse shape filter response, for the purpose? of 

step 40 of FIG. 3, is shown. IThe pulse shapesfilter response in the time- 
domain gi is inpufrinto an analog-to-digital conyeiyter sampling at a rate of 
{Ti*K2/Kih as shown in step 80^- .The digitised pulse sjiape filter response is 
25 then transformed into the frequency-domain by a fast fourier transform 

(FFT), in step 82, to produce the frequencyrdomain samples of the digitised 
pulse shape filter gi. The complex'gains gi (gi: i=l to K2) are then 
. ; normalised by breaking: the batch of K2 complex gains gi into K2/N batches 
-of?N;samples (gij: i=l to N, j=l . to K^/N); as -shown in, steps 84 to 88. 
30 The gains gij are squared and added to each other,. as in step 92, and 

an r arbitrary design parameter may be added to them to form the ; 
normalising parameters (qi:;i=l'to N) givei> by: qi=sum[ I gij l 2 +ei : j=l to 
K2/N, i=l to N], as in step 92.> The values for gij are normalised (divided) by 
the qi parameters as dn steps 94 to 98,<to generate -the K2 and ai 
35 parameters as shown in steps 100, 102 and 104. : . . 

Advantageously, a frequency-domain implementation is optimal for 
/ any received c pulse i.e. for pulses designed without the Nyquist ISI-free 
criterion or for pulses undergoing channel linear time-invariant filtering 
which introduces ISI. 
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The computational complexity of digital receiver operation is reduced 
by performing signal processing computations in the frequency-domain. The 
following method and explanation is, given by way of further explanation of 
the computational simplification achieved. 
5 The method describes the means by which the receiver can be 

designed and implemented in the frequency-domain using digital signal 
processing techniques, A frequency-domain implementation is especially 
suited for the TDMA case since an entire received slot can be analysed in 
"batch" processing while in conventional (continuous time) sequential 
10 processing is called for. 

A major advantage of the methodology is that it can provide optimal 
solutions with no' extra complexity to arbitrary received signal shape. As 
such it circum vents the fundamental Nyquist criterion for Inter Symbol 
Interference (ISI) free pulse shape (J.G. Proakis, Digital Communication, 
15 McGraw Hill). In the time-domain an equivalent receiver requires 

equalization in the form of a; (infinite-optimal or finite-sub-optimal) tap 
delay filter which entails heavy computational load. 

The implications of removing the Nyquist constraints are: 

(1) the pulse shape can be designed to accommodate other system 
20 considerations such as inter-channel interference (ICI), 

(2) the channel time invariant frequency distortion can be accommodated by 
recovering the symbols from the received pulse shape, absorbing the RF 
channel and IF filtering trivially into the receiver, and 

(3) filtering (e.g. to remove carrier-feedthroiigh) can be applied with little 
25 effect on performance and insignificant computational load. 

In the present method an alternative solution to equation [2] is 
presented. An apparatus is introduced in the frequency-domain which 
recovers the symbols from the received signal where ho special means (i.e. 
equalization) to overcome ISI is needed - it is inherent in the receiver 
30 structure. The basic equations derived below have their equivalent in the 
time-domain. However, for TDMA signals, analysis of an entire received 
slot using digital signal processing techniques enables the implementation 
of the receiver in the frequency-domain. ' '*' 

The mathematical expressions associated with the method are now 
35 presented, by way of completeness. 

The method describes the apph cation of a discrete fourier transform 
. of the received signal and the efficient estimation of the symbols 
{s(n):-0...N-l} from the transformed signal. ' 
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Applying a discrete fourier transform (BIT) to the received signal of 
equation [1], we obtain 



x((D) = g(co)X s ( n ) e 2T( _ j wT - i:i) ' 

5 . n ©eQ k ={© k .k = 2^k/K where : k = 0...K-l} [3] 

= g(cp)w^(cpT) T s\ . -. . . 

wlicrc . * - * * 

~w N (cpTf. = ^ . [4 ] 

.s^[s(0),s(l)--s(N-l)] , - .v 

where ' '-- 

10 co represents frequency, bold letters are used for variables in the frequency- 

' d6main 5 vectors are arrowed variables (e.g. s). and the -notation operator 

- . • {vector} implies the operation oneach element of the vector:. Furthermore, 

shorthand for vectors are used, as in equation {4]. 



15 Rewriting equation {3] in matrix form (matrices are in capital letters) 
x = G.W.s ■ • '■ •■■ 
where "•'* ■ .: . 

• x = [x(© k .0),x(cb k .l)r.^x(cOk.(K-l))f.. . . . 

diag[&<»i:oy i g<a> u sl),-..:iB(<»*.(K~ 1))].. (51a 

■ w = [w N (co k :o);, w N (co k . ij,.: : . w N (co k .(k - i))f ■ 

where diag, [ j is the diagonal matrix of the entries in the bracket. 
20 Ignoring Q[ ] to the symbol Alpha-Bet, the unconstrained optimal estimator 
(Least Squares and Maximum Likelihood estimator assuming white 
... Gaussian noise) is: . . ._. 

'/J. ''"I * iHG.w.r 1 * / ' * ''; _ [6] 

., where A; 1 t is the generalised (pseudo) inverse of A. 

For ejjuation ^}-tobe meaningful we have K >=2N, i.e. the number of 
points? in the fourier,transform of x(t) is greater, than twice the number of 
symbols, to ^e estimated,' Although this is not a must, it follows naturally 
from' the dimension of x(t) which must o>ey the Nyquist sampling rate. 
Thus equation [6] seemingly entails the inverse of a K by N matrix, a 
formidable task for a practical DSP base4* receiver. 

Assuming the received pulse shape is known, a computationally 
efficient solution of equation [6], utilising a FFT and an IFFT as a 
consequence of some (non-trivial) mathematical manipulations, is provided. 
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An aspect of the invention is the implementation of a computationally 
efficient solution of equation [6]. . 

If T is an integer, namely the sampling rate is a whole multiplicative 
factor of the symbol- rate and choosing K=NT in the DFT gives: 



[7] 



•W = [W N (1),W N (2)./.W N 
W N (1) = W N (2) =... W N = [^n (co N .0), w N (co N . 1),... w N (co N (N - 1))] 



Assuming K=NT is not a restriction, K should be chosen to be larger than 
NT so as not to truncate the first and last symbols' pulse shape. This is 
equivalent to augmenting the symbol vector by zero symbols on both frame 
sides. • . , ( . ' . ( . 

Applying the definition of the Generalised inverse to equation [6] and 
noting thafc-G is diagonal and that W is orthohormal we get: 

(GW) _1 x = (W H G H GW) _1 W H G H x 

.. ■ ^XW^PG.wJ.XwSGfx, 

vi=i ^ - - y i=i - : ...... 

. . =N.(W N ):^f G?G i j'SG-S, . ... 
• = NlIFFT I^^.GfGj.j 



where the sub-equations [enumerated] are explained as follows: 
sub-equation 0: definition. > I 

sub-equation 1: substitute equation [7] after breaking the matrices 
20 and vectors per the definition in equation [8]. 

sub-equation 2: obvious 
sub-equation 3: use orthogonal property * - 
sub-equation 4: definition. 



25 s= 1/N.IFFT 

where 



2\ 



-1 



( T 



V i = l 



[8] 
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Xj=i-th N dimensional sub vector of x, 
Gi=i-th N dimensional diagonal sub matrix of G, and 
is the hermitian transpose of G • 
' . • . ! - c ' ■ -■, ■■ 

Equation [8] is a key equation of the mathematical derivation. It 
implies that given the DFT of the data and the pulse shape G, estimating 
the symbols involves first multiplying, Gfii (K complex multiplications), 
secondly, overlapping and summing ]T G? • * i complex add- 

i=l 

multiplications) and thirdly applypg .a 1$ ppint IFFT. 



It. should be noted that 
(i) the external product by 



in equation [8] is not considered from 



a computational view point as it can be. absorbed as a normalisation factor 
into GH; 

• .- /it 1 ■ 2\7* 



15 (ii) for a Nyquist pulse we have 



^IGj] = a constant and thus 

.... A!?>- ! J . 

normalisation is not required (see J.G. Proakis, Digital Communication, 

McGrawHill); 

( x " 2 V 1 

(iii) if 



^|Gi| is ill conditioned (e.g. has large amplitude variation) the 

- i=1 J ■ ... 
symbol estimates will be nbise sensitive. 1 The remedy to this problem is well 

20 known and is reflected in the mean square error case, described later, 

(iv) the weighted least squares solution rather than least squares solution 

may be employed to accommodate noise variations across the spectrum, 

which results in modifying equation [8] to give 

if t -^V l -Y Tl " x " ' ; < M • '"• ' A 

Zl w Al..j::{•X w i G ^ i .)}• • 

25 where Wi is a real diagonal weighting matrix. . . 

The extension of the above to the case tha,t Tis a rational number is 
as follows. - 
Let 

30 T=Q/R - • ■ ' . ^ 



1 = 1/N.IFFT 
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where Q and R are integers. Then, if the signal is interpolated by a factor R 
the result of the previous section can be applied. However, interpolation by 
R in the frequency domain is trivially, achieved by zero padding the vector x 
by (R-l) K dimensional zero vectors (correctly placed on the unit circle). 
5 Applying equation [8] to the padded x is per the previous section. The 
actual padding is not needed as long as G is expressed in the correct 
interpolated sampling rate and : the overlap-add operation is done in the 
interpolated domain. 

Namely, the zero padding is not performed in practice, furthermore, 
10 note that 'the 3)FT of the sampled- data" is performed prior to the virtual 
interpolation explained above. 

The receiver expressed in eqiiation[8] and its extension for rational 
sampling rate in the symbol rate T is depictedliri FIG. 2. 

Using the mean square error criterion, rather than the least squares 
15 criterion used above, results in replacing equation [8] by : . 

f t • 2y r / T > * \ 

where N is the diagonal matrix of the noise spectral density vector. 

The frequency-domain receiver, just as the time-domain receiver is 
the optimal receiver in the least squares sense, as in equation [8], or mean 
20 square error sense, as in equation [10]. Since the time-domain and the 
frequency-domain representation are related by a linear operator, the 
. performance of the two receivers are identical for ISI-free signals. . 

A major benefit of the frequency-domain receiver is that it is not 
restricted to Nyquist pulses. Indeed, the solution in equations [8] and [10] 
25 are optimal for any received, pulse i.e., for pulses designed without the 

Nyquist ISI free .criterion and for pulses undergoing channel linear time- 
invariant filtering. Qr^ the other hand, ISI effects are usually compensated 
for in a time^domain receiver by, equalization - an expensive and not 
t necessarily optimal (if a finite, tap delay approach is used) operation. 
30 Moreover, having a complex channel impulse response requires a complex 
finite impulse response filter for equalization purposes which makes the 
solution in the time-domain even more expensive computationally. 

Thus an improved method of recovering symbols of a digitally 
modulated radio signstl, which avoids some of the disadvantages of prior art 
35 arrangements, is provided. . , 
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Claims 



1. A method of recovering a finite duration of N symbols of a digitally- 
modulated received radio signal comprising the steps of: 
5 . . receiving a signal having a symbol rate; 

sampling the signal at a sampling ;period Tj to generate a first batch 
of Ki time-domain samples* wherein 1/Ti is an integer multiple of the 
symbol fre.quency.of the received signal and Kl is chosen to be a multiple of 

N; , 

10 . transforming the batch of Ki time T domain samples to the frequency- 

domain thereby providing a first-batch of Kl frequency^domain samples; 

multiplying the first batch,of Ki fre.quency-domain.samples by a set of 
frequency-domain Kj. complex gains (ai: i=l :: to Kl), to rprovide a second 
batch of . Ki^equency^omain^samples;,, . 

15 breaking the second batch of Kl, frequency-domain samples into Ki/N , 

batches of N, samples; . ..• . 

adding the Ki/N batches of N samples to each other to generate a 

single batch of N samples and A - 

inverse transforming the single batch of N samples into the time- 
20 domain to obtain N 'samples (of desired symbols). 



2. A method according to' claim 1 i; wherein the Ki complex gains (ai: i=l 
^ to Kl ), are calculated by evaluating the received pulse' shape filter response 

' in the freduency-domain to obtain the batch of Kl complex frequency- 
25 domain samples ! (gi: i=l to' K'l), wherein the parameters gi represent the 
' frequency-domain pulse shape filter response. 

3/ a method of recovering a fu^te duration of N symbols of a digitally- 

r modulated- received radio signal comprising the steps of: 

30 "*' receiving a i signal having a symbol rate; 

sampling the iignal at a sampling period Ti to generate a batch of Kl 

' time^doinain samples, wherein 1/Ti is not an integer multiple of the symbol 
frequency of the received signal; 

' transforming the batch of Ki time-domain samples to the frequency- 
35 domain thereby providinga batch of Ki frequency-domain samples; 

appending zeros to the batch of Ki frequency-domain samples, 
thereby providing a batch of K2 samples in the frequency-domain (wherein 
K2 > Ki) , such that K2= (T2 *N) samples, 1 wherein LT2 is an integer 
multiple of the symbol frequency; and 
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multiplying the batch of K2 frequency-domain samples by a set of 
frequency-domain K2 complex gains (ai: i=l to K2), to provide a batch of K2 
frequency-domain samples. 

5 4. A method according to claim 3, wherein the K2 complex gains (ai: i=l 
to K2), are calculated by evaluating the received pulse shape filter response 
in the frequency-domain to obtain the batch of K2 complex frequency- 
domain samples (gi: i-1 to K2), wherein the parameters gi represent the 
frequency-domain pulse shape filter response. v 

10 

5. A method according to any of the preceding claims wherein the 
parameters gi representing the frequency-domain pulse shape filter 
response are pre-stored. 

15 6. A method according to claims 3 and 4 wherein the complex gains gi 

(gi: i=l to K2) are normalised, comprising the "steps of: 
\ breaking the batci of K2 complex gains gi into K2/N batches of N 

samples (gij: i=l to N,j=l to K2/N); : . 

adding the absolute squares of the K2/N batches to each other to 
20 generate a single batch of N normalising values (qi: i=l to N) such that 
qi=sum[lgij|2 : j =lto K2/N];- 

normalising (gij: i=l to N) for each j by dividing the values gij by the 
values qi to obtain a new set of gains (aij=aij/qi : i=l to N, j=l to K2/N); 

combining the K2/N subgroups (aij) back into a single group (ai) to be 
25 used as the complex gains of ai in claim 1. 

7. A method according to claim 6 wherein to each qi there is added a 
constant ei as a design parameter. 

30 
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